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Abstract 

We approximately compute the normalization constant of the first infrared renor- 
malon of the pole mass (and the singlet static potential). Estimates of higher order 
terms in the perturbative relation between the pole mass and the MS mass (and in 
the relation between the singlet static potential and a s ) are given. We define a match- 
ing scheme (the renormalon subtracted scheme) between QCD and any effective field 
theory with heavy quarks where, besides the usual perturbative matching, the first 
renormalon in the Borel plane of the pole mass is subtracted. A determination of the 
bottom MS quark mass from the T(IS') system is performed with this new scheme and 
the errors studied. Our result reads m bMs( m &Ms) = 4 210lgg (theory)!^ (a s ) MeV. 
Using the mass difference between the B and D meson, we also obtain a value for 
the charm quark mass: ™ c Ms( m c Ms) = ^ ^ICT^q (theory)^!? (m b ms) +45(^1) MeV. We 
finally discuss upon eventual improvements of these determinations. 

PACS numbers: 14.65.Fy, 14.65.Dw, 12.38.Cy, 12.39.Hg 



pineda@part icle . uni-karlsruhe . de 



1 Introduction 



Systems composed by heavy quarks are very important in the study of the QCD dynamics. 
This is due to the fact that they can test QCD in a kinematical regime otherwise unreachable 
with only light quarks. They are characterized by having an scale, the mass of the heavy 
quark, much larger than any other dynamical scale in the problem. Therefore, it seems 
reasonable to study these systems by using effective field theories where the mass has been 
used as an expansion parameter. Some examples of such are HQET IJ for the the one-heavy 
quark sector and NRQCD § or pNRQCD § g for the Q-Q sector. 

On general grounds, the matching coefficients of the (QCD) effective theories suffer from 
renormalon ambiguities |J (for a review on renormalons see ||). This means that, in prin- 
ciple, one can not compute these matching coefficients with infinity accuracy in terms of 
the (short distance physics) parameters of the underlying theory. From a formal point of 
view, there is not fundamental problem related with these ambiguities, since the renormalon 
ambiguities of the matching coefficients cancel with the renormalon ambiguities in the cal- 
culation of the matrix elements in the effective theory in such a way that observables are 
renormalon free (as they should be) up to the order in the expansion parameters to which 
the calculation has been performed. 

Being more specific, a generic matching coefficient c[yjm) would have the following 
perturbative expansion in a s : 



oo 

C ( Z y/ m )=C+^C rl < +1 , (1) 
n=0 



where a s (in the MS scheme) is normalized at the scale v. Its Borel transform would be 

oo ±n 

iW) = £c n - (2) 
and c is written in terms of its Borel transform as 

oo 

c = c + J dte~ t/a ° B[c}{t). (3) 
o 

The ambiguity in the matching coefficient reflects in polesf] in the Borel transform. If we 
take the one closest to the origin, 

5B[c](t) ~ (4) 

where a is a positive number, it sets up the maximal accuracy with which one can obtain 
the matching coefficients from a perturbative calculation, which is (roughly) of the order of 

8c~r n *af, (5) 

where n* ~ — . Moreover, the fact that a is positive means that, even after Borel resumma- 
tion, c suffers from a non-perturbative ambiguity of order 

5c ~ (Aqcd)^ • (6) 
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Observable 2 



Figure 1: Symbolic relation between observables through the determination of the matching 
coefficients of the effective field theory. 

As we have mentioned, these renormalon ambiguities will cancel anyhow in the final cal- 
culation of the observable. Therefore, we could ask ourselves why bother about this problem. 
The answer comes from the procedure we use to relate (and then predict) observables (see 
Fig. [I]). Schematically, one takes one observable to fix the matching coefficient. In order to 
relate the matching coefficient with the observable, we use the effective field theory as a tool 
to provide the power counting rules in the calculation. Since we are relating one observable 
(a renormalon free object) with a matching coefficient suffering of some renormalon ambi- 
guity, there must be another source of renormalon ambiguity as to cancel this one. As we 
have mentioned, the latter comes from the calculation of the matrix element in the effective 
theory. There is a problem here, however. If the matrix element in the effective theory is a 
non-perturbative object, the ambiguity of the matching coefficient is of the same size than 
this non-perturbative object that we do not know how to calculate anyhow. Nevertheless, 
if the matrix element in the effective theory is a perturbative object, the only way it has 
to show the renormalon is in a bad perturbative behavior in the expansion parameters in 
the effective theory (what is happening is that the coefficients that multiply the expansion 
parameters are not of 0(1) due to the renormalon, i.e. the renormalon is breaking down 
the assumption of naturalness implicit in any effective theory). This may seem irrelevant 
since the ambiguity is the same in either case but in the latter situation it means that the 
observable is less sensitive to long distance than the matching coefficient itself. Thus, if, for 
instance, we wanted later on to get the short distance parameters from that (weakly sensitive 
to long distance physics) observable, we are not doing an optimal job, since we use an inter- 
mediate parameter (the matching coefficient) that can not be obtained with better accuracy 
than the ones displayed in Eqs. whereas the observable (and the short distance 

parameter) is less sensitive to long distance physics (the same problem will also appear if we 
want to relate two weakly sensitive to long distance physics observables through a more long 
distance sensitive matching coefficient). We can consider three examples (observables) to 
illustrate this point: the mass of the meson B (Mb), the inclusive semileptonic decay width 
of the B (T(B — > Xuluj) and the mass of T(IS') (assuming in this latter case ma 2 3> Aqcd 
for illustration). For these observables, the first non-perturbative corrections (leaving aside 
renormalons) are of the following type: 



where mos is the pole mass. The above results only become true if the perturbative piece 
can be computed with such precision. Nevertheless, this is not true in the on-shell (OS) 

1 In general, this pole becomes a branch point singularity but this does not affect the argumentation. 



A 4 



QCD 



5 np M B ~ Aqcd 



S np T(B -> XJu) ~ G 



' 2 177 3 \ 2 




(7) 
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scheme, where one uses the pole mass as an expansion parameter, since the latter suffers 
from renormalon ambiguities [^, §. Therefore, effectively, the above observables can only be 
computed perturbatively (working with the pole mass) with the following precision 

5%^M B ~ A QCD , 5^T(B -> XJv) ~ G|m 4 os A QCD , <^ ert -)M T(ls) ~ Aqcd- 

(8) 

For the first observable this discussion is irrelevant but not for the other two, which become 
rather less well known that they could be. 

Roughly speaking, the above discussion means that the renormalon of the matching 
coefficients can be spurious (it is not related to a real non-perturbative contribution in the 
observable) or real (it is related to a real non-perturbative contribution in the observable). 
In fact, this distinction depends on the observable we are considering rather than on the 
renormalon of the matching coefficient itself^. The point we want to stress is that at the 
matching calculation level it makes no sense this distinction. Therefore, there is no necessity 
to keep the renormalon ambiguity (that it only appears due to the specific factorization 
prescription we are using) in the matching coefficients (even more if we take into account 
that its only role is to worsen the perturbative expansion in the matching coefficients). Thus, 
our proposal is that one should figure out a matching scheme where the renormalon ambiguity 
is subtracted from the matching coefficients. This is the program we will pursue here for 
the specific case of effective field theories with heavy quarks. In principle, the same program 
should, eventually, be carried out with other effective theories. Nevertheless, the worsening 
of the perturbative expansion was especially evident in effective theories with heavy quarks. 
This is due to the fact that the renormalon singularities lie close together to the origin and 
that perturbative calculations have gone very far in this case P|-||22|. In other effective field 



theories this problem may have not become so acute (yet) but it may become relevant in the 
future. 

The structure of the paper is as follows. In the next section, we compute the normalization 
constant of the first infrared (IR) renormalon of the pole mass. We also give estimates of the 
higher order coefficients of the perturbative series relating the pole mass with the MS mass. 
In section [|, we compute the normalization constant of the first IR renormalon of the singlet 
static potential and also estimates of the higher order perturbative terms in the potential 
are given. In section f|, new definitions of the pole mass and the singlet static potential are 
given, within an effective field theory perspective, by subtracting their closest singularities 
to the origen in the Borel plane. In section [5], we provide a determination of the bottom 
MS mass from the T(15) mass and an estimate of the errors within this new approach. In 
section || we provide a determination of the charm MS mass by using the mass difference 
between the B and D mesons. In section [7], we give our conclusions and discuss how to 
improve our determinations of the bottom and charm masses. 



2 In the discussion above the renormalon of the pole mass would be called real if we think of the B meson 
mass and spurious if we think of the mass of the T(1S). 



3 



2 Mass normalization constant 



The pole mass can be related to the MS renormalized mass - which in principle can be 
measured to any accuracy at a very high energy scale - by the series 

oo 

m os = ™ms + J2 r ««" +1 > ( 9 ) 

n=0 

where the normalization point v = ttlj^ is understood for (in this way we effectively 
resum logs that are not associated to the renormalon since both m^si^) and m Ms( m Ms) 

do 

not suffer from the bad renormalon behavior) and the first three coefficients tq, r\ and r 2 are 
known || (a s = a^ l \ where ni is the number of light fermions; we will assume in this work 
that any other quark except the heavy quark is massless). The pole mass is also known to 
be IR finite and scheme- independent at any finite order in a s |23] . We then define the Borel 
transform 

°? oo j.n 

mos = mm + J dte-^ B[m os }{t), B[m os }{t) = £ r n — (fO) 

The behavior of the perturbative expansion of Eq. at large orders is dictated by the closest 
singularity to the origin of its Borel transform, which happens to be located at t = 27r//3 , 
where we define 

Being more precise, the behavior of the Borel transform near the closest singularity at the 

M 

4-7T 



Iv- 



or igin reads (we define u — — " 



B[mos)(t(u)) = N m v ^_\ u y +i (l + c i(! - 2 «) + c 2 (l - 2u) 2 + ••■) + (analytic term), 

(11) 

where by analytic term, we mean a piece that we expect it to be analytic up to the next 
renormalon (u — 1). This dictates the behavior of the perturbative expansion at large orders 
to be 

™„ (PoY r(rc+l + i>) / 6 b(b-l) \ 

r " = N ""[ii) m+b) ( 1 + (^n? Cl + (n + t) (n . + t-i) C2 + ---J- (12) 

The different b, ci, c 2 , etc ... can be obtained from the procedure used in [p4j| . The coefficients 
b and c\ were computed (exactly!) in p3{] , and c 2 in (where, apparently, there are some 
missprints for its analogous expression s 2 ). They read 



2$ 



b = wk, (13) 
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and 

° 2 " 6(6 = 1) 32^| • (15) 

We then use the idea of [EH] (see also [EH]) and define the new function 



D m (u) = Y. D t ) u n = {l-2u) 1+b B[m os \{t{u)) (16) 

n=0 

= N m u (l + ci(l - 2u) + c 2 (l - 2m) 2 H ) + (1 - 2m) 1+b (analytic term) . 

This function is singular but bounded at the first IR renormalon. Therefore, we can expect 
to obtain an approximate determination of N m if we know the first coefficients of the series 
in u and by using 

N m v = D m {u= 1/2). (17) 

The first three coefficients: D$, D$ and Dff are known in our case. In order the calculation 
to make sense, we choose v ~ m. In other words, we avoid to have another large (small) 
parameter, otherwise we should dealt with the necessity of resummation. If, for illustration, 
we restrict ourselves to the large (3$ approximation, the dimensionful parameters rearrange 
in the quantity 

/ v \ 2u v 
m — ~ v{\ + (2u - 1) In 1 }. 

\mj m 

Therefore, the underlying assumption is that we are in a regime where (besides 2u — 1 <S 1) 

(2m- l)ln— < 1. 
m 

For the specific choice v = m, we obtain (up to 0(u 3 )\ u= i/ 2 ) 

N m = 0.424413 + 0.137858 + 0.0127029 = 0.574974 (n f = 3) (18) 
= 0.424413 + 0.127505 + 0.000360952 = 0.552279 (n f = 4) 
= 0.424413 + 0.119930 - 0.0207998 = 0.523543 (n f = 5) 

The convergence is surprisingly good. One can also see that the scale dependence is quite 

mild (there even appears a place of minimal sensitive to the scale dependence, see Fig. ||). 
If, for illustration, we take — — = 2 one obtains, for nf = 4, N m = 0.545580. On the other 

m MS 

hand, for small — — , the scale dependence starts to become important (if, for illustration, 

m MS 

we take ^ = 1/2, one obtains, for n f = 4, N m = 0.434619). 



By using Eq. (|T2|), we can now go backwards and give some estimates for the r n . They are 
displayed in Table |l| We can see that they go closer to the exact values of r n when increasing 
n. This makes us feel confident that we are near the asymptotic regime dominated by the 
first IR renormalon and that for higher n our predictions will become an accurate estimate 
of the exact values. In fact, they are quite compatible with the results obtained by other 
methods like the large (3q approximation (see Table [I]). If we also compare with the estimates 



for r3 given in Ref . [28 , our results are also roughly compatible with those up to differences 



of the order of 10 %, 20 % and 30 % for n/ = 3, n/ = 4 and nj = 5, respectively. 
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Figure 2: x = dependence of N m for nf = 4. 



In our estimates in Table [TJ, we have included (formally) subleading terms in the 1/n 
expansion up to 0(l/n 2 ). We show their impact in Table 0. They happen to be corrections 
with respect the leading order result in all cases and of the order of 4 % or smaller for n > 2. 
The smallness of the 0(l/n) corrections with respect the leading order result can be traced 
back to the approximate (numerical?) pattern f3 n ~ (3q (for (3\ and ^2) for nj = and that 
nf terms do not jeopardize this rule (this explains the strong dependence of the results with 
nf). In fact, if this rule were exact, one can easily see that c n = for all n. The breaking 
of this rule explains the non-zero values for c n . The 0(l/n 2 ) terms depend on /?3, which has 
been recently evaluated in [29|, through C2. In this case, however, we rather have f3^ ~ 2(3q 
and the nj terms are quite large breaking the pattern above. In fact, c\ (the upper entry 
for r in Table 0) strongly depends on nf so that if for Uf = 0, C\ and Ci (the lower entry 
for r in Table §) are of the same size [c\ ~ —0.215, c<i ~ 0.185), for nf = 5, C\ becomes 
an order of magnitude smaller than c<i- Therefore, one may think of the smallness of C\ as 
a numerical accident (for some specific values of nj) not reflecting the natural size of the 
0(l/n) terms and therefore explaining the apparent breakdown (or slow convergence) of the 
1/n expansion. This would be elucidated if ^4 were known. In the mean time, we will stick 
to this belief. 

We can now try to see how the large (3 approximation works in the determination of N m . 
In order to do so, we study the one chain approximation from which we obtain the value M 



jy(large /3 ) 



c 



f £ 



0.976564. 



7T 



(19) 



By comparing with Eq. (|18|), we can see that it does not provide an accurate determination 
of N m . This may seem to be in contradiction with the accurate values that the large (3q 
approximation provides for the r n (starting at n = 2) in Table [I]. Lacking of any physical 
explanation for this fact, it may just be considered to be a numerical accident. In fact, the 
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Tn — r n/ m MS 


^0 


h 


r-i 




r 4 


exact (rif = 3) 
Eq. (0) (n f = 3) 
large (3 (n f = 3) 


0.424413 
0.617148 
0.424413 


1.04556 
0.977493 
1.42442 


3.75086 
3.76832 
3.83641 


18.6697 
17.1286 


118.441 
97.5872 


exact (rif — 4) 
Eq. (0) (n, = 4) 
large (3 (n/ = 4) 


0.424413 
0.645181 
0.424413 


0.940051 
0.848362 
1.31891 


3.03854 
3.03913 
3.28911 


13.8151 
13.5972 


80.5776 
71.7295 


exact (n/ = 5) 
Eq. (|D ( n/ = 5) 
large (3 (n f = 5) 


0.424413 
0.706913 
0.424413 


0.834538 
0.713994 
1.21339 


2.36832 
2.36440 
2.78390 


9.73117 
10.5880 


51.5952 
51.3865 



Table 1: Values of r n for v = mj^g. Either the exact result (when available), the estimate 
using Eq. or the estimate using the large (3q approximation 



— r n/ m MS 


^0 


fi 


r 2 


r 3 


f 4 


0(l/n) (n f = 3) 
0(l/n 2 ) (n f = 3) 


-0.164 


-0.046 


-0.027 


-0.019 


-0.015 


0.237 


-0.103 


-0.017 


-0.007 


-0.004 


0(l/n) (n f = 4) 


-0.105 


-0.028 


-0.016 


-0.012 


-0.009 


0(l/n 2 ) ( n/ = 4) 


0.274 


-0.126 


-0.020 


-0.008 


-0.004 


0(l/n) (n f = 5) 


0.024 


0.006 


0.003 


0.002 


0.002 


0(l/n 2 ) (n f = 5) 


0.326 


-0.165 


-0.023 


-0.009 


-0.005 



Table 2: 0(1 /n) corrections (normalized with respect the leading solution) of our r n estimates 
for different number of light fermions. 



agreement between our determination and the large (3q results does not hold at very high 
orders in the perturbative expansion, whereas we believe, on physical grounds, since our 
approach incorporates the exact nature of the renormalon, that our determination should go 
closer to the exact result at high orders in perturbation theory. Nevertheless, the large /3q 
approximation remains accurate up to relative high orders. 



3 Static singlet potential normalization constant 

One can think of playing the same game with the singlet static potential in the situation 
where A QCD C 1/r. The potential, however, is not an IR safe object at any order in the 
perturbative expansion |30], |TTJ. Its perturbative expansion reads 

oo 

K (0) (r;^) = EK {0 n ) < +1 , (20) 

n=0 

where we have made explicit its dependence in the IR cutoff v us . The first three coefficients 
V}q, v}x and v}°2 are known JT(| as well as the leading- log terms of v} ^ |ll| (for the 
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renormalization-group improved expression see fl2]|). Nevertheless, these leading logs are 
not associated to the first IR renormalon, since they also appear in momentum space (see 
also the discussion below), so we will not consider them further in this section (or if one 
prefers, we will only consider coefficients that we fully know). 

We now use the observation that the first IR renormalon of the singlet static potential 
cancels with the renormalon of (twice) the pole mass. This has been proven in the (one-chain) 
large @ approximation in |3T], |j] and at any loop (disregarding eventual effects due to u us ) in 
f32fl . We would like to argue that this cancellation holds without resort to any diagrammatic 
analysis as far as factorization between the different scales in the physical system is achieved. 
This can be done within an effective field theory framework where any renormalon ambiguity 
should cancel between operators and matching coefficients. Therefore, in the situation where 
1/r ^> Aqcd, one can do the matching between NRQCD and pNRQCD and 2mos + V s ^ 
can be understood as an observable up to 0(r 2 Aq CD , AQ CD /m) renormalon (and/or non- 
perturbative) contributions (see Eq. (0)) 3 . This would prove the (first IR) renormalon 
cancellation at any loop (as well as proving the independence of this IR renormalon of 
v us ). In a way, the argumentation would be similar to the one made for the renormalon 
cancellation in the electromagnetic correlator using the operator product expansion |?|, 
where the renormalons are absorbed in local condensates. In our case, however, we are 
talking of non-local objects, but fortunately, effective field theories provide themselves as 
useful in these situations. 

We can now read the asymptotic behavior of the static potential from the one of the pole 
mass and work analogously to the previous section. We define the Borel transform 



K (0) 



^ oo ±n 

/ die-'/- B[V™](t) , B[V®](t) = £ ygf- (21) 



The closest singularity to the origen is located at t — 2^/(3$. This dictates the behavior of 
the perturbative expansion at large orders to be 

v mn^ ((3 \ n T(n + l + b) ( b 6(6-1) \ 

V « n ~ NvU {^) T(l + b) l v 1+ (^T6) Cl+ (n + 6)(n + 6-l) C2 + "j' (22) 

and the Borel transform near the singularity reads 

B[V s (0) }{t(u)) = N v u- ]— — (l + ci(l - 2u) + c 2 (l - 2uf + •••) + (analytic term). 



;i - 2u) 



(23) 

In this case, by analytic term, we mean an analytic function up to the next IR renormalon 
at u = 3/2 |3§. 

As in the previous section, we define the new function 

oo 

D v (u) = ^^"Hl-^^K^iM) (24) 

n=0 



3 Note that the same argumentation does not apply for the octet static potential V , . The reason is that, 
even at leading order in 1/m, 2mos + is not an observable. This is due to the fact that there is still 
interaction with low energy gluons, as one can see from Eq. (Q). Therefore, one expects 2mos + Vo°^ to be 
ambiguous by an amount of O(Aqcd)- 



8 



= N v u (l + ci(l - 2u) + c 2 (l - 2m) 2 + •••) + (1 - 2m) 1+6 (analytic term) 

and try to obtain an approximate determination of Ny by using the first three (known) 
coefficients of this series. By a discussion analogous to the one in the previous section, we 
fix v — 1/r. We obtain (up to 0(u 3 ) 



\u=l/2) 

Ny = -1.33333 + 0.571943 -0.345222 = -1.10661 
= -1.33333 + 0.585401 - 0.329356 = -1.07729 
= -1.33333 + 0.586817 - 0.295238 = -1.04175 



3) 
4) 
5). 



(25) 



The convergence is not as good as in the previous section. Nevertheless, it is quite acceptable 
and, in this case, apparently, we have a sign alternating series. In fact, the scale dependence is 
quite mild (see Fig. ||]) except (again) for small values of vr. Overall, up to small differences, 
the same picture than for N m applies. 



•1 . 05 
-1 . 1 
•1 . 15 
-1.2 
•1 . 25 
-1.3 
1 . 35 
-1 . 4 



N 



v 



0.60.8 1 1.21.4 1.61.8 2 

x 



Figure 3: x = vr dependence of Ny for nj = 4. 

So far we have not made use of the fact that 2N m + Ny = 0. We use this equality as a 
check of the reliability of our calculation. We can see that the cancellation is quite dramatic. 
We obtain 



2N m + N v 



2N m - N 



v 



0.038 
0.025 
0.005 



n f = 3 
n f =A 

n f = 5. 



We should stress that the evaluation of N m and Ny uses independent inputs. Therefore, 
this cancellation seems to be nontrivial. This makes us confident that the number we have 
obtained for N m (or Ny) is quite accurate. The difference is even better than what one 
would have expected from the last terms in the series of Ny, but in this case, this may due 
to the fact that the series is sign-alternating. 
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In the following we will use the determination of N m to fix Ny. Any difference should be 
included in the errors. 

We can now obtain estimates for Vj® by using Eq. fl2"2"|). They are displayed in Table 
|3|. Note that in Table |3| no input from the static potential has been used since even Ny 
have been fixed by using the equality 2N m = —Ny. We can see that the exact results are 
reproduced fairly well (the same discussion than for the r n determination applies). This 
makes us feel confident that we are near the asymptotic regime dominated by the first IR 
renormalon and that for higher n our predictions will become an accurate prediction of the 
exact results. The comparison with the values obtained with the large (3q approximation 
would go (roughly) along the same lines than for the mass case, although the large (3q results 
seem to be less accurate in this case (see Table |3|) . 



T/(0) _ y(0) 
* s,n * s,n 


v s,0 


v s,l 




v® 




exact (n/ = 3) 
Eq. (|D {n f = 3) 
large (3 (n f = 3) 


-1.33333 
-1.23430 
-1.33333 


-1.84512 
-1.95499 
-2.69395 


-7.28304 
-7.53665 
-7.69303 


-37.3395 
-34.0562 


-236.882 


exact (jif = 4) 
Eq. © (n/ = 4) 
large (3 (n f = 4) 


-1.33333 
-1.29036 
-1.33333 


-1.64557 
-1.69672 
-2.49440 


-5.94978 
-6.07826 
-6.59553 


-27.6301 
-27.0349 


-161.155 


exact (nf = 5) 
Eq. (P ( n/ = 5) 
large /3 (n f = 5) 


-1.33333 
-1.41383 
-1.33333 


-1.44602 
-1.42799 
-2.29485 


-4.70095 
-4.72881 
-5.58246 


-19.4623 
-21.0518 


-103.190 



Table 3: Values ofVjp^ with v = 1/r. Either the exact result (when available), the estimate 
using Eq. $22), or the estimate using the large (5q approximation ^Tl\ , |77 



In order to avoid large corrections from terms depending on u us , the predictions should 
be understood with v us = 1/r and later on one can use the renormalization group equations 
for the static potential to keep track of the v us dependence. 

Finally, we would also like to mention that the previous discussion about the determina- 
tion of N m in the large (3q approximation also applies here for the determination of Ny. 



4 Renormalon subtracted matching and power count- 
ing 

In effective theories with heavy quarks, the inverse of the heavy quark mass becomes one of 
the expansion parameters (and matching coefficients). A natural choice in the past (within 
the infinitely many possible definitions of the mass) has been the pole mass because it is 
the natural definition in OS processes where the particles finally measured in the detectors 
correspond to the fields in the Lagrangian (as in QED). Unfortunately, this is not the case in 
QCD and one reflection of this fact is that the pole mass suffers from renormalon singularities. 
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Moreover, these renormalon singularities lie close together to the origin and perturbative 
calculations have gone very far for systems with heavy quarks. At the practical level, this 
has reflected in the worsening of the perturbative expansion in processes where the pole 
mass was used as an expansion parameter ]13|, |2(|. It is then natural to try to define a 
new expansion parameter replacing the pole mass but still being an adequate definition for 
threshold problems. This idea is not new and has already been pursued in the literature, 
where several definitions have arisen. For instance, the kinetic mass ||, the PS mass [32 1 



the IS mass and the PS mass |[35| . We can not resist the tentation of trying our own 
definition. We believe that, having a different systematics than the other definitions, it 
could further help to estimate the errors in the more recent determinations of the MS quark 
mass. Our definition, as the definitions above, try to cancel the bad perturbative behavior 
associated to the renormalon. On the other hand, we would like to understand this problem 
within an effective field theory perspective. From this point of view what one is seeing is 
that the coefficients multiplying the (small) expansion parameters in the effective theory 
calculation are not of natural size (of 0(1)). The natural answer to this problem is that we 
are not properly separating scales in our effective theory and some effects from small scales 
are incorporated in the matching coefficients. These small scales are dynamically generated 
in n-loop calculations (n being large) and are of 0(me~ n ) (we are having in mind a large (3q 
evaluation) producing the bad (renormalon associated) perturbative behavior. The natural 
way to deal with this problem would be to perform an expansion of the (small) scale m e~ n 
over the (large) scale m, or alike, for a n-loop calculation. Unfortunately, this is something 
that standard dimensional regularization does not know how to achieve since it does not 
know how to separate the scale m from the scale m e~ n treating them on the same footing^. 
In order to overcome this problem, we may think of doing the Borel transform. In that case, 
the renormalon singularities correspond to the non-analytic terms in 1 — 2m. These terms also 
exist in the effective theory. Therefore, our procedure will be to subtract the pure renormalon 
contribution in the new mass definition, which we will call renormalon subtracted (RS) mass, 
mRS (with no pretentious aims, all the other mass definitions do cancel the renormalon as 
well, but rather for notational purposes). We define the Borel transform of mRs as follows 



B[m RS } = B[m os ] - N m v f -——^ b (l + c x {l - 2u) + c 2 (l - 2uf + •••), (26) 



4 This problem does not appear in hard cut-off renormalization schemes or alike (to which, after all, will 
belong our scheme), which explicitely cut-off these scales. At this respect, we can not avoid thinking that 
we are in a similar situation to the beginnings of NRQCD (NRQED), where it was better known how to 
achieve the separation of scales with hard cut-off than with dimensional regularization |3(| . This problem 
was first solved in (3?], || within an effective field theory framework (see also |3£| where a solution within 
a diagrammatic approach was provided). A key point in the solution came by realizing that the way to 
implement the separation of scales (matching) in dimensional regularization was by first expanding the 
Feynman integrals with respect the small scales (the ones that would be kept in the effective theory) prior to 
integration. Obviously, a similar solution here (with no necessity of hard cutoff), where one should expand 
with respect the small scale me~ n prior to integration, would be most welcome. This would change the 
perturbative expansion and the (small) scale me~ n should be kept explicit in the effective theory. To date, 
we are not able to further substantiate this discussion but we expect to come back to this issue in the future. 
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where Uf could be understood as a factorization scale between QCD and NRQCD and, at 
this stage, should be smaller than m. The expression for reads 

( \ X-M (P°Y n+lf T(n + l + b-k) 

where Cq — 1. We expect that with this renormalon free definition the coefficients multiplying 
the expansion parameters in the effective theory calculation will have a natural size and also 
the coefficients multiplying the powers of a s in the perturbative expansion relating with 
m^fg. Therefore, we do not loose accuracy if we first obtain m^s and later on we use the 
perturbative relation between m^s and in order to obtain the latter. Nevertheless, since 
we will work order by order in a s in the relation between m^g and m^, it is important to 
expand everything in terms of a s , in particular a s (z/y), in order to achieve the renormalon 
cancellation order by order in a s . Then, the perturbative expansion in terms of the MS mass 
reads 

oo 

™RsO/) = ™Ms + E (28) 

n=0 

where r^ s = (my^, v,Vf). These r^ s are the ones expected to be of natural size (or at 
least not to be artificially enlarged by the first IR renormalon). 

From Eq. (P7|), we can see that we are subtracting a finite piece (that should correspond 
to the renormalon for large n) for every r n . It is more than debatable whether we should 
do any subtraction for tq. Therefore, to test the scheme dependence of our results, we also 
define a modified RS scheme, RS', as follows 

( \ X-N (&>Y n + x ( ^ T(n + l + b-k) 

mRsivf) = m os - 2^N m v f a s + (u f ) ^ c k ^ - - - ^ , (29) 

and the Borel transform corresponds to 

B[m RS ,) = B[m RS ] + N m v f (1 + c x + c 2 + ■ ■ •) . (30) 

At this stage, we would like to make some preliminary numbers to estimate the effect of 
our definition as well as to compare with other threshold masses. In order to simplify the 
discussion we define an static version of the IS mass: 

m (stat 1C ) ^ + V^r) = m _ + / _ ^ + _ _ _ ^ ^ 

and compare the expansions for the pole mass, the PS mass, the RS mass and the (static) 
IS mass. We display the results in Table |] for the bottom quark and in Table |5| for the top 
quark. By default the value a s (M z ) = 0.118 is understood. The four loop evolution equation 
has been used for the running of a s as provided by the program RunDec.m [RS]. Overall, 



all the threshold mass definitions work well in the cancellation of the renormalon, as we can 
see by comparing with the pole mass result at higher orders. For the bottom quark case, we 
have displayed results for Uf = 2 and z// = 1 since on the one hand, on conceptual grounds, 
we would like to keep Uf below (or of the order of) what it will (roughly) be the typical 
scales of the inverse Bohr radius in the T(IS') system but, on the other hand, we would like 
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Masses 


0{a s ) 




OK 3 ) 


o{<4) 


total 




401 


199 


144 


147 


5102 


™rs 

MRS' 

m PS 

( static ) 


111 
401 
210 
102 


50 
114 
80 
50 


17 
38 
42 
19 


7 
15 

8 


4 395 
4.778 
4 542 
4 389 


fn RS 
m R& 
m PS 

(static) 


256 
401 
306 
251 


95 
157 
120 
94 


40 
74 
67 
41 


21 
41 

22 


4 622 
4.882 
4.703 
4619 



Table 4: Contributions at various orders in a s for different mass definitions for the bottom 
quark case, either with Uf = 1/r = 2 GeV (middle panel) or with Uf — 1/r = 1 GeV (lower 
panel). The results are displayed in MeV. For the O(af) results, the estimate from Table \^ 
has been used. The other parameters have been fixed to the values ^ms(^ms) = 4.21 GeV, 

V = m Ms( m Ms) and n f = A - 



Masses 


0{a s ) 


0(a%) 


0{al) 


OK) 


total 




7.585 


1.615 


0.497 


0.221 


174.917 


™RS 

«2 RS ' 

m PS 

(static) 

m\ s 


7.355 
7.585 
7.447 
7.368 


1.469 
1.590 
1.518 
1.475 


0.391 
0.461 
0.460 
0.395 


0.139 
0.180 

0.141 


174.354 
174.816 
174.425 
174.379 



Table 5: Contributions at various orders in a s for different mass definitions for the top 
quark case. The results are displayed in GeV. For the 0{a\) results, the estimate from 
Table [| has been used. We have used the following set of parameters: Uf — 1/r — 3 GeV, 
m Ms( m Ms) = 165 GeV,v = m I g(m S g) andn f = 5. 

to keep it significantly larger than Aqcd- This gives us little room to play. For the IS (static 
mass), we have chosen, for a closer comparison with the other definitions, 1/r = 2, 1 GeV, 
respectively. For the top quark case, we have chosen Uf — 1/r — 3 GeV. 

The shift from the pole mass to the RS mass affects the explicit expression of the effective 
Lagrangians. In particular, in HQET, at leading order, a residual mass term appears in the 
Lagrangian 

C = h (iD - 5m RS ) h + O (—) , (32) 

where Sm^s — ^os — m RS and similarly for the NRQCD Lagrangian. 

For heavy quark-ant iquark systems in the situation where Aqcd "C ma s , it is convenient 
to integrate out the soft scale (~ ma s ) in NRQCD ending up in pNRQCD. If we consider 
the leading order in 1/m, the residual mass term is absorbed in the static potential (in going 
from NRQCD to pNRQCD, one runs down the scale Uf up to uj < ma s ). We can then, 
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analogously to the RS mass, define an singlet static RS (RS') potential 



VSsCRff) M = V r + 25m RS( Rso , (33) 

where the coefficients multiplying the perturbative series should be of 0(1) (provided that we 
expand t^ ' and Sm^s in the same parameter, namely a s ). Notice also the trivial fact that 
the scheme dependence of m RS cancels with the scheme dependence of Vrs- It is interesting 
to see the impact of this definition in the improvement of the perturbative expansion in the 
potential. Then, following analogously the discussion for the masses, we have compared our 
definition (the RS potential) with the PS potential and the singlet static potential for some 
typical values appearing in the bottom and top quark case. We have displayed the results in 
Tables || and |7|. We see that the improvement is quite dramatic with respect the expansion in 
the singlet static potential, yet, again, for the bottom case, we have little room for changing 



Potentials 


0(a s ) 


OK) 


o(o») 


o«) 


total 


K (o) 


-910 


-306 


-302 


-383 


-1902 


v s,RS 


-205 


3 


-2 


-3 


-208 


v s,RS' 
V s,PS 


-910 
-446 


-54 
-42 


-14 
-25 


-6 


-984 
-513 


y(0) 

y(0) 
V s,RS' 

v s,PS 


-558 
-910 
-678 


-63 
-180 
-116 


-41 

-95 
-75 


-26 
-54 


-687 
-1239 
-869 



Table 6: Contributions at various orders in a s for different singlet static potential definitions 
for some typical scales in the T system, either with Vj = 2 GeV (middle panel) or with Vf = 1 
GeV (lower panel). The results are displayed in MeV. For the 0{oQ) results, the estimate 
from Tabled has been used. The other parameters have been fixed to the values v = 1/r = 2.5 
GeV and nf = 4. 



Potentials 


0(a s ) 


0(ai) 


O(ai) 




total 


V W 


-5 679 


-874 


-404 


-237 


-7194 


y(0) 
v s,RS 

v s,RS' 
y(0) 


-5 077 
-5 679 
-5 317 


-548 
-788 
-648 


-185 
-294 
-233 


-67 
-120 


-5 876 
-6 881 
-6199 



Table 7: Contributions at various orders in a s for different singlet static potential definitions 
for some typical scales in t-i systems near threshold. The results are displayed in MeV. For 
the O(a^) results, the estimate from Tabled has been used. The parameters have been fixed 
to the values v = 1/r = 30 GeV, Vf = 3 GeV and rif = 5 . 
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The pNRQCD Lagrangian in the RS scheme formally reads equal than in the OS scheme: 



j pNRQCD — J d 3 Hd 3 r 



'M«'l«b-^- + ^l~^(0-^-^ + ...)s 



+Ot ( iDo -^--V%(r) + ..)o 
\ m RS J j 

Vb{t), 



+gV A (r)Tv {O f r • E S + S f r • E O} + £-y^Tr {O f r • E O + O f Or • e} 

d 3 R~G%G^ a , (34) 

where V}][ s = Vj°) +25mRs(RS') and the 1/m potentials also get (straightforwardly) affected 
by rewriting the expansion in 1/mos hi terms of I/uirs (see |3|, |4], [L4[] for the definitions in 
the OS scheme and details). The above Lagrangian provides the appropriate description of 
systems for which, for their typical r, one has the inequality 1/r ^> Aqcd- If one further 
assumes that r 2 AQ CD <C ma 2 ., the power counting rules tell us that the leading solution 
corresponds to a Coulomb-type bound state being the non-perturbative effects corrections 
with respect to this leading solution. We will assume to be in this situation in the following. 

One of the ultraviolet cutoffs in pNRQCD is u us , which cutoffs the three-momentum of 
the gluons in pNRQCD and fulfills the relation mot 2 s <C v us <C ma s . Without any further 
assumption about the (perturbative or non-perturbative) behavior at scales below u us , one 
can compute the heavy quarkonium spectrum. Formally, we have the following expression 
(at the practical level, one could work in the OS scheme and do the replacement to the RS 
scheme, with the proper power counting, at the end) 

oo 

M nlj = 2m RS + £ 4jf S K>r + SM^(u us ) , (35) 

m=2 

where the v us scale dependence of the different pieces cancels in the overall sum (for the 
perturbative sum this dependence first appears in ). 

We expect that by working with the RS scheme (or with any other achieving the renor- 
malon cancellation) the coefficients multiplying the powers of a s will now be of natural size 
and therefore the convergence improved compared with the OS scheme. So far, the A™'^ s 
coefficients are exactly known for n = 2, 3, 4 whereas some partial information is also known 
for A^ s and 5M^ S \ We will discuss them further in the next section. 

At this stage, we would like to discuss some few theoretical issues that may appear in the 
readers mind with respect these investigations (see also the discussion in Ref. P0| ). First, 
once one agrees to give up using the pole mass as an expansion parameter, one may still 
wonder why not to use the MS mass instead. There are several answers to this question. If 
we consider pNRQCD, working with m^g would mean introducing a large shiftQ, of 0(ma s ), 
in the pNRQCD Lagrangian, therefore jeopardizing the power counting rules. In fact, the 

5 This is certainly so for t-i physics. Nevertheless, for the bottom, the 0(ma s ) term does not seem to 
be that large numerically (see Table ^), being much smaller than the typical values of the soft scale in the 
T(15 f ). Therefore, it may happen that working with the MS mass does not destroy the power counting 
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same happens in HQET. Thus, it seems necessary to work first in an scheme where the 
power counting rules are preserved, this invalidates the TO^g mass (as far as we are talking 
about processes where the heavy quarks are near OS). On the other hand, we also require 
not to have large (not natural) coefficients multiplying the perturbative expansion. This 
invalidates the pole mass, remaining only the threshold schemes. Nevertheless, once one 
observable is obtained in a (safe) threshold scheme, one could consider to rewrite it in terms 
of the MS mass. Even that could be eventually dangerous. If, for instance, we consider the 
quantity M n ij and expand it in terms of to^, there appear, in principle, two problems. One 
has to do with the resummation of logs, which, after all, is one of the motivations of the 
whole factorization program between different scales that effective field theories are (another 
is to provide, in an easy way, the power counting rules of the dynamics). By expanding 
everything in terms of a s , we reintroduce a potentially large log, hi 22 ), in the coefficients 
multiplying the powers of a s (note that we can not minimize this log if at the price of 
introducing another large log, ln 2222 ^). Another problem is that, due to the fact that there 
is another scale, ma s , besides to, at least on conceptual grounds, we would not achieve the 
renormalon cancellation order by order in a s but it would occur between different orders in 
a s jeopardizing, in principle, the convergence of the perturbative expansion. 

From the above discussion, we have seen that it is crucial to have a renormalon free 
expansion parameter and an scheme that preserves the power counting rules or, equivalently, 
to be able to write an effective Lagrangian (with the dynamics, i.e. the power counting rules) 
in terms of the new threshold masses. This is something that it is achieved by the threshold 
schemes, i.e. the kinetic, the PS-like, the IS and the RS schemes. However, the IS scheme 
seems to rely on assuming that the Q-Q system is a mainly perturbative system. Note also 
that the IS and PS schemes depend on u us . We see that the power counting rules in the 
RS scheme are the same than the ones in the OS scheme, the difference being that it is 
now expected that the terms in the perturbative expansion will be of natural size. We have 
then been able to solve the renormalon problem without giving up the factorization between 
different scales that is provided by effective field theories (in the OS scheme) . Note also that 
the RS mass only knows of scales above the cut-off of the effective theory, a desirable feature 
in any factorization program. 

Throughout this work much emphasis has been put on working with effective field theo- 
ries. Therefore, one may honestly ask whether the problems with renormalons here exposed 
would disappear if one gives up using effective field theories. Nevertheless, a closer inspection 
seems to show that the renormalon problem always appears in physical problems where one 
has different scales and wants to achieve factorization between them. Thus, any procedure 
(name it effective field theory or not) that provides factorization and power counting rules 
will have the same kind of problems. 



rules of pNRQCD (or HQET) at the practical level. However, it docs from a conceptual point of view, in 
particular, one would have problems to assign power counting rules. Moreover, it would not incorporate 
the expected physical fact that for scales between m and i>f » Aqcd the threshold masses are the proper 
expansion parameters for processes near on-shell. 
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5 Bottom MS quark mass determination 



In this section, we will determine the MS (and RS) bottom mass from the T(IS') mass. We 



will use the known results at next-to-next-to-leading order (NNLO) in the OS scheme |13| and 
rewrite them in the RS scheme where the coefficients multiplying the perturbative expansion 
are expected to be of a natural size (for the moment we neglect ultrasoft contributions). In 
Figures ^ and |5|, we plot the scale dependence of the LO, NLO and NNLO predictions for 
the T(1S) mass in the RS and RS' scheme. 
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Figure 4: We plot 2mfcRg (dashed line), and the LO (dot-dashed line), NLO (dotted line) 
and NNLO (solid line) predictions for the T(1S) mass in terms of v in the RS scheme. The 
value of m^Rs is taken from Eq. ffidj ). 

We can also check the dependence of our results with respect to the theoretical and 
experimental parameters. As a first estimation of the errors, we allow for a variation of v, 
u f , a s and N m as follows: v = 2.5^- 5 GeV, v f = 2±l GeV, a s (M z ) = 0.118 ± 0.003 and 
N m = 0.552 ± 0.0552. For the RS scheme, we obtain the following result^ 



m 6iRS (2 GeV) = 4 387j| 8 



i6V<Xs)+lt(N m ) MeV; 



m^im,^) = 4 203+2 5 (i/);5( Z / / )- 2s 



v — 10/ 



(36) 

"6,ms \"%,M8J ~ ^^ wlJ +25v;+6V/;+27V"s; + io(^V m ) MeV. (37) 
For the RS' scheme, we obtain the result (with the same variation of the parameters) 

m 6 , RS v(2 GeV) = 4 782;^(z/); 3 7 (z/ / )l^(« s );l(V m ) MeV; (38) 

m fe,Ms( m b,Ms) = 42 1jgW(»//)3(tt i a(jy w ) MeV. (39) 

6 Here and in the following, in the determination of ni^jg, we have used our estimate of the four-loop 
relation. 
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Figure 5: We plot 2m feRS / (dashed line), and the LO (dot-dashed line), NLO (dotted line) 
and NNLO (solid line) predictions for the T(1S) mass in terms of v in the RS' scheme. The 
value of m^RS' is taken from Eq. fiSdj). 



The MS mass result depends weakly on the variation of Vf and N m . For the latter it may 
seem surprising in view of the strong dependence of the RS mass but, to some extent, the 
variation of N m can be understood as a change of scheme. 

With the results obtained in Eqs. (|36|) and (0), the expansion for the Y(IS') mass would 
read in the RS scheme (see Eq. (j35[)): 



M T(15) = 8 774 + 559 + 120 + 7 MeV, 



and in the RS' scheme: 



T(15) 



9 564- 158 + 56 - 2 MeV 



(40) 



(41) 



Both series seem to show convergence. Nevertheless, only in the RS' scheme a more physical 
interpretation can be given where the leading order solution corresponds to the negative 
Coulomb binding energy (but we can see that this affects very little the determination of the 
m^jg mass) . 

Naively, from Eqs. (|4"0"D-(4"I), one would be very optimistic about the errors. Nevertheless, 
let us try to go deeper into the error analysis. First, it may be not realistic to conclude from 
these results that the magnitude of the NNLO terms is of order of few MeV. The relative 
size between different orders will depend on the scale (as one can see in Figs. [| and [5]). The 
major problems come at small v and we will discuss them later. Leaving them aside, for the 
RS' scheme, the NNLO contributions remain quite small, whereas for RS scheme they can 
be of the order of —50 MeV for large values of v. Being even more conservative, one can 
split the NNLO contributions into the ones due to the singlet static potential (+62 MeV and 
45 MeV, with v = 2.5 GeV, for the RS and RS' scheme, respectively) and the relativistic 
(—55 MeV and —47 MeV, with v = 2.5 GeV, for the RS and RS' scheme, respectively) 
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ones. Each of them happen to be of the order of 50 MeV with opposite signs. As far as 
the relativistic corrections is concerned, we can see that they are comfortably smaller than 
the 0{ma 2 s ) leading term. For the correction due to the singlet static potential, convergence 
is also found in the RS scheme, whereas in the RS' scheme the convergence is quite slow. 
In fact, this is due to the fact that the magnitude of the LO and NLO terms in the RS 
scheme is larger than in the RS' scheme. Nevertheless, the discussion on the magnitude of 
the corrections due to the static potential is quite delicate since it strongly depends on the 
value of Vf used. A closer look also shows that the main part of the scale dependence comes 
from the relativistic corrections, which become very important at small scales. Nevertheless, 
one can think of this scale dependence arising because of neglecting some higher order logs 
that, on theoretical grounds, we know anyway. The, 0(mosctg), NNNLO logs are known in 
the OS scheme |TJJ (see also |TR [HI]). They may be considered to be of two different origen. 



Either those which can be minimized by choosing v of the order of the inverse Bohr radius 
(for these an explicit expression can be found in Ref. ^^), or those related either to the hard 
(~ mm), or to the ultrasoft scale (~ lnz/ us ). Therefore, the only piece of Al^ s left unknown 
(although by far the most difficult one) is the log-independent one. For this, there exists, 
at least, a large /3 determination [|Tj], In order, this evaluation in the OS scheme to be 
useful, one has to change to the RS scheme, where the renormalon cancellation is explicit 
and the NNNLO correction is expected to be of natural size (note that one has to take into 
account the change from the pole mass to the RS mass in the expressions). Once a value 
for the log-independent term of A^ s is assumed (we will use the large (3q result), one can 
compute Al^ s without any further ambiguity and, thus, the complete NNNLO correction. 
We will use this result to give an estimate of the NNNLO corrections (and, therefore, further 
substantiate our previous discussion on the convergence of the series), and in order to study 
the scale dependence. Note that at this stage a dependence on v us appears. If we believe 
that the large /?o approximation provides (numerically) a good estimate of, at least, the size 
of the renormalon contribution to the binding energy in the OS scheme, we could expect our 
determination to provide a rough estimate of higher order effects. We show our results in 
Figs. | and % 

Let us first concentrate on the scale dependence. We first consider the result without 
the inclusion of the logs related to the hard/ultrasoft scale (being them of different physical 
origen). The results are the dot-dashed line in Figs. ^| and [7| At small scales, we see 
the typical oscillatory behavior between different orders in the perturbative expansion when 
the series ends to converge. For scales where the calculation is reliable, we see that the 
correction goes in the direction one would expect by choosing a somewhat smaller value for 
v (and therefore closer to the scale mRsC/a s ~ 2 GeV) in the NNLO evaluation. On the 
other hand, if we are concerned about the absolute size of the NNNLO corrections, what we 
can see is that, according to our estimate, they are small. We can now add the large logs 
associated to the hard and ultrasoft scale (we use Eq. (19) from Ref. [JT3] ). They depend 
on v us . We have chosen v us = 1 GeV. Our results correspond to the dot lines in Figs. |6] 
and 0. We see that the correction is relatively small (and it goes closer to the dot-dashed 
result when we increase v us ). It also blows up for small values of v (this makes us believe 
that in order to properly deal with the scale dependence at small v a renormalization-group 
improved evaluation could be useful). Therefore, we feel relatively confident about these 
sources of NNNLO corrections. However, this evaluation can not give an estimate of genuine 
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0(mRs«s) corrections coming from other sources. We will discuss some of them right now 
as well as non-perturbative effects. 
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Figure 6: We 'plot the NNLO (dashed line) versus the NNNLO estimates in the RS scheme, 
a) NNNLO estimate without hard/ultras oft logs (dot-dashed line) and b) NNNLO estimate 
including hard/ultras oft logs with v us = 1 GeV (dotted line). The value of taken 
from Eq. (^). 



So far, we have just considered the expansion in A™'^ s a™ and neglected the ultrasoft 
effects. Let us consider them now. An explicit expression for the leading correction due 
to ultrasoft effects can be obtained using the multipole expansion without any assumption 
about the perturbative or non-perturbative behavior at ultrasoft scales. It reads |14| 



(we actually write the Euclidean expression for easier comparison with other results in the 
literature) 



5M„ 



us, 

nl 1 



Tf_ 

3N n 



where Hj} s 



■ 



™RS 



^(n,/|re- i ^ RS -^ S )r|n,/)^E a (t)0(t,O):^E b (O))(z/ us ), 
o 

(42) 

— mRsC/«s/(4n 2 ). Different possibilities appear 



and E^ s 



2N C r 

depending on the relative size of Aqcd with respect to the ultrasoft scale ma 2 s . Let us first 
study whether the T(IS') lives in the situation where Aqcd "C ma 2 s . In this case, the scale 
ma 2 can be treated perturbatively and one can lowers further the cutoff up to the situation 



Aqcd < ^ 



s <C ma 2 . Scales below u' us can be parameterized in terms of local condensates 

A 2 



within an expansion in O ( Ja^S ) • Eq. (fE|) would then read (with leading- log accuracy) 



5M. 



US/ 
nl \ 



SM VS, pert. + SM m, no-pert 



nl 



(43) 



5M US,pert. 



E^ln 

n m RS a s 



^ + 4C A Cj 



n(2l + l) 



2C 2 



n(2l + V 



n- 
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Figure 7: We plot the NNLO (dashed line) versus the NNNLO estimates in the RS' scheme, 
a) NNNLO estimate without hard/ultras oft logs (dot-dashed line) and b) NNNLO estimate 
including hard/ultras oft logs with v us = 1 GeV (dotted line). The value of m&^s' is taken 
from Eq. (jgfty. 



3n 



and 



6M VS, no-pert. 



c n T: 



(44) 



n=0 



1/ (m 



RS 



) and O n ~ Aqqtj (see [^2| for details). CqOq was first computed 



and CiOi in Ref. 



The perturbative leading-log v us scale dependence would 



where C, 
in Refs. 

cancel against the scale dependence of Al'^ s producing log(a s )-type terms. The validity of 
these results for the physical system under study rely on two assumptions that one should 
check. One the one hand, the final result will depend on a s (ma 2 s ) (this is easily seen within 
an effective field theory framework) but ma 2 s is an small scale for the T(1S). If we take as 
an estimate (H — Ei) w ~ 360 MeV, we see that perturbation theory is not truthsworthy for 
this scale (although, obviously, numerical factors can play a role). If one, anyway, wants to 
study this situation, one can include the log-dependent terms from Eq. ( ^3|) in the NNNLO 
perturbative estimate (the dependence on v us effectively disappears). This would produce 
a ~ +50 MeV shift in the RS scheme (~ +59 MeV in the RS' scheme) with respect to the 
NNLO results (for v = 2.5 GeV), becoming even larger for small scales but decreasing for 
larger values of v (in any case the shift depends on finite piece prescriptions) with, somewhat, 
the same shape than the dotted lines in Figs. |6| and 0. On the other hand, as far as the non- 
perturbative corrections is concerned, one should check that the operator product expansion 
in terms of local condensates converges. If we use our central values (in the RS' scheme) 



plus the central values for the condensates used in Refs. [13], |42 |, we obtain 



C O + C x O x = 144 - 143 MeV . 



(45) 
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We see that we do not find convergence. The situation improves by lowering the scale v in a s 
(97 - 66 MeV for v = 2 GeV and 53 - 21 MeV for v = 1.5 GeV) and it also depends on the 
values of the condensates, which are poorly knownQ We see that this working hypothesis 
is not favored by the central set of parameters of the T(IS') system, although it can not 
be ruled out if one scan over the possible values of the parameters (in fact, if one considers 
v — 2 GeV a more natural scale for the soft scale, the operator product expansion would be 
on the verge of convergence). Nevertheless, for t-i production near threshold, the situation 
Aqcd "C ma 2 may be applicable and the results explained above useful. 

It may then seem that the T(IS') system lives in the situation where Aqcd ~ ma 2 s . In 
that case, we can not lower further v us and the chromoelectric correlator in Eq. ( fi2]) cannot 
be computed using perturbation theory. We are then faced with the necessity of computing 
a non-perturbative non-local condensate. If in the previous case the knowledge of the local 
condensates was poor the situation is even worse now. This non-local condensate is related 



to the gluonic correlator for which the most general parameterization reads []44 

(gF^ixMz, 0)^gF b pa (0)) = (5^ - 5M [v(x 2 ) + V,(x 2 

dVi(x 2 ) 

~f~ (%v%<jfi[ip •Ev%p3[i<j •Efi%(r3vp 4" •Efi'Epfii/a) ~Qg2 ' ^ ^) 

where for the gauge string a straight line is understood and D and T>\ are invariant functions 
of x 2 . In our case the following combination appears^] 

(sE a (t)0(t, 0)^E 6 (0)> = 3 (v(t 2 ) + V.it 2 ) + t 2 ^p-) . (49) 



In Ref. [ £7|] , a lattice evaluation of the gluonic correlator was performed. The following 
parameterization was used to describe the lattice data (with t > 0) 

V(t 2 ) = ^ exp (-*/A ) + A exp {-t/\ A ) V 1 {t 2 ) = ^ exp (-*/A ) + M exp {-t/\ A ) . 

(50) 



7 The difference with the conciusions in Ref. [|42], Q follows from the fact that in these works a s , as 
defined in fl3|| , was used as the expansion parameter instead of a s used here. Unfortunately, a s suffers from 
the renormalon ambiguity making it potentially large and bringing convergence to the operator product 
expansion. 

8 At this stage, we would like to report some discrepancies with the results of |45|| . If we consider the 
situation mv 3> Aqcd S> niv 2 , the leading non-perturbative effects can be parameterized in terms of a 
potential term as follows jl6| |J 



rp fOO 

SV = r 2 ^ J o dt( ff E«(t)0(i,O):^E h (O))K s ). 



(47) 



In Ref. Ea] , the following result was reported: 



SV = ^ rdt(v(t 2 ) + lvi(t 2 )). (48) 



N C J V v ' 2 



We have difficulties to accommodate this result with Eq. p7| ) after using Eq. (pf9[). In particular, Eq. ( ^q ) 
does not appear to be able to reproduce the leading logs predicted by perturbation theory. 
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One would expect the 1/i 4 terms to have something to do with the perturbative terms. In 
any case, at short distances, the gluonic correlator behavior should go closer to the one 
expected by perturbation theory. Unfortunately, we see no indication of this but rather 
the opposite. Whereas the perturbative result predicts a behavior ~ — 1/t 4 for the (E(t) • 
E(0)) operator, the lattice simulations get a positive slope at short distance (either for 
quenched or unquenched simulations). Lacking of any explanation for this fact, we will 
refrain of using those results in order to get an estimate of the non-perturbative effects 
(another point of concern is that the gluon condensate prediction for quenched simulations 
is one order of magnitude larger than the phenomenological value). Therefore, we will not 
try to give any number for the ultrasoft contribution in this paper and add them to the 
errors^. In order not to overestimate them, we will constraint the allowed range of values 
for the ultrasoft contribution by some consistency arguments. If we rely on power counting 
rules, since Aq C d ~ ma 1i the non-perturbative effects would formally be of NNLO and from 
the discussion above about the perturbative NNLO effects one would assume a value around 
~ 50 MeV. On the other hand, for consistency of the theory (AQ CD /(ma s ) 2 <C ma 2 ), the 
ultrasoft corrections should be smaller than the leading order solution (of order ~ 160 MeV). 
Note that if this consistency argument is not fulfilled, the very same assumption that we can 
describe, in first approximation, the T(IS') by a Coulomb-type bound state fails. In order to 
keep ourselves as conservative as possible, we would only demand the ultrasoft corrections 
to be smaller than the binding energy. Therefore, we will assign to our evaluation of the 
RS mass a ±50 MeV error. This is roughly equivalent to assign a ±100 MeV error to the 
evaluation of the binding energy from the ultrasoft contributions. This error may also be 
considered compatible with the discussion of the situation where ma 2 3> Aq C d- From the 
previous discussion about the magnitude of the NNNLO contributions, we will add another 
±25 MeV error to our evaluation of the RS mass (roughly equivalent to assign a ±50 MeV 
error to the evaluation of the binding energy) from any other (perturbative) source of higher 
order effects. This may seem conservative in view of the estimates of the NNNLO effects 
shown in Figs. ^| and [7| and the related discussion but we also include the expected error 
from finite charm mass effects (see |18|). Therefore, from the RS scheme evaluation, our 
final estimate for the bottom mass reads (in order to avoid double counting with the error 
from higher order effects, we do not include now the scale dependence error) 

m b , RS (2 GeV) = 4 387^(US + higher orders)^(a a );g(JV m + u f ) MeV; (51) 

™ 6j msK,ms) = 4203j$(US + higher orders) +%(a,)+%(N m + u f ) MeV. (52) 

Whereas from the RS' scheme evaluation, we obtain the result (with the same variation of 
the parameters) 

m 6iRS '(2 GeV) = 4 782+£j(US + higher orders) t^(a 8 )^(N m + v f ) MeV; (53) 

™&,msK,ms) = 4 214j$(US + higher orders) +MM+\l{N m + u f ) MeV. (54) 
We average the two values obtained for the MS mass. We then obtain (rounding) 

%msK,ms) - 4210^(theory);iK) MeV, (55) 

9 In principle, a similar problem may appear in sum rules calculation. See the discussion at the end of 
this section. 
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where we have added (conservatively) the theoretical errors linearly. These figures compare 
favorably (within errors) with other determinations of the bottom mass. Either with other 
determinations using the T(1S) mass (the second references in [ 2"T| , [52]] and ref. flIB|), sum 
rules | 2~T| , [18 1, lattice or from measurements at LEP [3j| . 



In this paper, we have just taken into account the first IR renormalon of the pole mass and 
the singlet static potential. Nevertheless, there are subdominant renormalons that eventually 
could play a role. On the singlet static potential side, one expects the first problems to come 
from a 0(AQ CD r 2 ) IR renormalon. For the pole mass there are 0(AQ CD /m) renormalons. A 
priori, it is not clear which one will play the dominant role. In fact, it will depend on the 
relative size between Aqcd and ma 2 s . 

In the situation where a description in terms of local condensates is appropriated [ma 2 ^> 
Aqcd), the leading genuine non-perturbative corrections to the mass scale like ?7?,(Aqcd / 'ma s ) A 
but this quantity is (parametrically) much smaller than the non-perturbative effects asso- 
ciated to the subleading renormalons, either from the pole mass or from the singlet static 
potential. Moreover, in this case, parametrically, the leading ambiguity would come from 
the subleading pole mass renormalon of 0(AQ CD /m). Therefore, we would be in a similar 
situation that when working with the pole mass, where the actual accuracy of the result 
was set by the perturbative calculation. Thus, one should first get rid of these subleading 
renormalons in order to improve the accuracy of the calculation before doing any reference 
about genuine non-perturbative effects. 

In the situation where ma 2 ~ Aqcd, the subleading renormalon ambiguities from the 
singlet static potential and the pole mass are parametrically of the same order than the 
genuine non-perturbative corrections. The Y(IS') seems to live closer to this situation. 

It is usually claimed that the non-perturbative effects in sum rules are smaller than in 
the Y(IS') mass. We would like to mention that, at least parametrically, this is not the 
case under the standard counting ~ a s (where n labels the moments in sum rules, see 
PH for details). Nevertheless, it may happen that they are numerically suppressed. This 
is indeed the case if one considers that one can describe the non-perturbative effects by local 
condensates [Q. However, one can only use the expression in terms of local condensates 
if one is in the situation where — ^> Aqcd- This would be analogous to the assumption 
ma 2 . ^> Aqcd and it may difficult to fulfill. Therefore, it is more likely that the non- 
perturbative corrections will also depend on a nonlocal condensate, in fact, arising from 
the same effective theory, on the same chromoelectric correlator that the T(IS') mass does. 
Thus, in order to estimate the non-perturbative errors in sum rules evaluations, it would 
be most welcome to have, at least, the explicit expression of the non-perturbative effects 
in the situation — ~ Aqcd, which, by now, is lacking. In that way, one could relate the 
non-perturbative effects for different moments in the sum rules or with the non-perturbative 
effects of the T(IS') mass and eventually search for less sensitive non-perturbative physics 
observables. 
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6 Charm MS quark mass determination 



In this section, we give a determination of the MS (and RS) charm mass. In principle, one 
could think of using the same procedure than in the previous section for the J/\I/ (or i] c ). 
This would imply to believe that the J/\I/ is a mainly perturbative system. We prefer to 
avoid this assumption in this work and to obtain the charm mass in a different way. Our, 
maybe weaker, assumption will be to assume that HQET can be used for charm physics. 
We will search for observables that are weakly sensitive to non-perturbative physics. The 
observable that we will choose is the difference between the B and D meson mass. Here, 
it follows an analogous discussion than in the introduction. Whereas both the B meson 
mass and the D meson mass are O(Aqcd) sensitive, the difference is 0(AQ CD /m) sensitive. 
Nevertheless, this improvement is not exploited if one writes this observable in terms of the 
pole masses. Therefore, here, we will use the RS masses and our previous determination of 
the bottom mass in order to obtain the charm mass. 

The (spin-averaged) mass difference between the B and D meson mass reads 

(M B ) - (M D ) = m fejRS - m c , RS + ^ ( — ^ — !— ] + 0(l/m RS ) , (56) 

where Ai can be related with the expectation value of the kinetic energy in HQET and 

(m b(d)) = Mb ^ ±^m , (57 ) 



The value of Ai is poorly known. We take the value Ai = 0.3 ± 0.2 (see [|5T| and references 
therein). We can now obtain the value of the charm mass. We perform the evaluation in 
both schemes, the RS and RS'. In order to estimate the errors, we fix Vf = 1 GeV and allow 
for a variation of a s , N m , Ai and nib as follows: a s (M z ) = 0.118±0.003, N m = 0.552±0.0552, 
Ai = 0.3 ±0.2 and m bME = 4 203^ MeV (in the RS scheme) and m bM g = 4 214±|g MeV (in 
the RS' scheme). For the variation of a s and N m , we use the correlated values of the bottom 
mass (that is why we do not include these errors in the variation of the bottom mass). For 
the error due to z/j see the discussion below. We obtain 

m CjRS (l GeV) = 1 181^(m biM g)^(« s );^(^ m ); 6 1(A 1 ) MeV; (58) 

MsKms) = 1 206« 6 7 K M s)^(a s )+^(iV m );l(A 1 ) MeV, (59) 
and for the RS' scheme 

m CiRS '(l GeV) = 1 477tS(m 6iH5 )^}(a.);g(JV m )3(Ai) MeV; (60) 

MsKms) = 1 207^K jM s)^(a s )^^(V m );i(A 1 ) MeV. (61) 

We can check that the perturbative relation between the RS and MS charm mass is indeed 
convergent. We obtain 

m CjRS (l GeV) = 1 206 - 53 + 20 + 6 + 3 = 1 181 MeV, (62) 
m C)RS '(lGeV) = 1207 + 205 + 46 + 13 + 6 = 1477 MeV. 
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In order to test our evaluation, since, for the charm quark, Uf stays close to the charm 
mass value (maybe jeopardizing the real structure of the perturbative expansion), we will also 
perform the calculation in the OS scheme (note that in that case large logs may appear). We 
use Eq. fl56[ ) with the replacement RS — > OS. In order to achieve the renormalon cancellation, 
we rewrite both m^os and m Cj os in terms of the respective MS masses and expand them in 
a s {m c jyfg). In this way, we obtain the renormalon cancellation for m^os ~ ^c,os but not for 
the 1/m terms. Therefore, for the latter, we use the MS masses as the expansion parameters. 
This effectively increases the magnitude of the 1/m terms but the other consistent option, 
to use the OS masses, is heavily affected by the renormalon (for instance m C) os — 2 263 MeV 
for m c jyjg = 1 210 MeV) introducing even larger errors and not reflecting the real size of the 
1/m corrections. We allow for the same variation of a s and Ai than above, whereas for the 
bottom mass we use m b y^ = 4 210 ± 90 MeV. We obtain 

«\msK,ms) = 1 254^K,Ms)il2(« s );^(Ai) MeV . (63) 

We can check the convergence of the perturbative expansion of m^os — m c,os i n terms of 
the MS masses in order to test the validity of this evaluation. We obtain (order by order in 

m btOS - m CiOS = 2 956 + 490 - 14 - 32 + 22 = 3 423 MeV . (64) 

For the last terms in the series the situation is not conclusive. One could think that the 
expansion is reliable up to, maybe, a ~ ±30 MeV uncertainty. One source of error comes 
from the expansion in a s since two kind of logs arise: In (ra&/ v) and In (m c / u), which can not 
be minimized at the same time. In fact, for the choice v = 2m c ^, we obtain m C Ms( m C Ms) = 
1 239 MeV and the expansion seems to improve: 

m bt0 s - m CiOS = 2 971 + 345 + 79 + 19 + 10 = 3 424 MeV. 

Let us consider further sources of error in the RS scheme evaluations. For these, we are 
using the RS bottom masses at a quite low Vf = 1 GeV. This produces that the convergence 
in the conversion from the MS to the RS and RS' masses for the bottom quark is slower, in 
particular for the RS' mass (see Table |]). One can then believe that higher orders in the 
relation between the MS and the RS' masses will add further positive contribution that in 
turn will increase the value of the MS charm mass bringing it closer to the MS evaluation. 
In any case, they are compatible within errors. Therefore, we add a ±40(20) MeV error 
to our RS'(RS) evaluation from the conversion from the MS to the RS'(RS) bottom quark 
mass. Another source of error comes from the 1/m terms. We have two effects compensating 
each other. With the RS scheme, the magnitude of the corrections in the relation between 
the RS and the MS bottom quark mass is smaller but at the price of worsening the 1/m 
expansion (this observation could also apply to the OS calculation). With the RS' scheme, 
the magnitude of the corrections in the relation between the RS and MS bottom quark mass 
is larger but the 1/m expansion improves (as well as being less sensitive to Ai). We estimate 
the 0(l/m 2 ) corrections to be of order ~ 15(30) MeV for the RS'(RS) evaluation and also 
add them to the errors. Overall, the theoretical errors in the RS and RS' evaluations are 
basically equivalent. Therefore, we, somewhat, weigh more the RS' scheme evaluation, being 
less sensitive to Ai. Our final result for the MS charm mass reads (rounding) 

™ c ,msKms) = 1 210^(theory)^K )Mg );«(A 1 ) MeV, (65) 
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where we have not included the errors in a s being negligible compared with the other sources 
of error. This result can be compared with two recent evaluations where Charmonium data 
was used 



At this stage, we can also give a prediction for A by using 

Ars = (M B ) - m biRS - + 0(1 /m 2 b RS ) . (66) 

We obtain (using m 6 ^g = 4 210 MeV) 

A RS (1 GeV) = 659 MeV, A M '(1 GeV) = 401 MeV. (67) 

We can see that it is crucial to specify the scheme in order to give a meaningful prediction 
for A. 

Finally, we would like to mention that m^Rs — m c,RS ( an d rn^os — m Cj os) suffers from 
renormalon ambiguities of 0(AQ CD /m) that cancel with the renormalon of Ai/m. Therefore, 
one could argue whether it makes any sense to give a value of Ai without specifying how to 
handle the 0(Aq CD /m) renormalon in m^Rg — m c ^. This could only be explained if the 
ambiguity due to the renormalon is much smaller than the genuine non-perturbative effects. 
If we take as an indication the perturbative expansions found above, we may believe that 
any ambiguity in the perturbative expansion is smaller than the genuine non-perturbative 
effects. On the other hand, this ambiguity may explain the spread of values one can find in 
the literature for Ai. 



7 Conclusions and outlook 

We have approximately computed the normalization constant of the first infrared renormalon 
of the pole mass (and the singlet static potential). Estimates of the higher order coefficients 
of the perturbative series relating the pole mass with the MS mass (and the singlet static 
potential with a s ) have been obtained without relying on the large (3q approximation. New, 
renormalon free, definitions of the mass and potential have been given, within an effective 
field theory perspective, by subtracting their closest singularities in the Borel plane. We 
have obtained the bottom MS quark mass from the T(IS') mass, and an estimate of the 
errors, within this new scheme. We have also obtained the charm MS mass by using the 
mass difference between the B and D mesons. Our final figures read 

™ 6j msKms) = 4 210^(theory);^(« s ) MeV (68) 

and 

™ c ,ms(™ c ,ms) = 1 210±™(theory)^K^);^(A 1 ) MeV. (69) 
Several lines of research may follow from our results. 

The use of conformal mapping could eventually lead to an improved convergent series in 
the evaluations of N m and Ny. 
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It would be interesting to apply the RS scheme to bottomonium sum rules or to t-t 
production near threshold and see how large the differences are with respect other determi- 
nations available in the literature. 

Our result for the bottom mass has been obtained in the zero mass charm approximation. 
Finite charm mass effects (IB], |5B| should be incorporated in future studies. 



We have seen that our evaluation suffers from scale dependence for small v. At NNLO, the 
main source of scale dependence comes from the relativistic corrections. The incorporation 
of NNNLO effects does not seem to correct this fact. It may happen that a renormalization- 
group improved result could solve this problem. It is worth noting that there already exists 
one result in the OS scheme in the situation where ma 2 s 3> Aqcd [0- It would be desirable 
to have an independent evaluation within pNRQCD and without relying on the unequality 
ma 2 3> Aqcd- After that, one should transform the renormalization-group improved results 
from the OS to a renormalon free scheme. This, a priori, may not turn out to be completely 
trivial. 

One of the (potentially) major source of errors in our evaluation of the bottom mass is the 
non-perturbative contribution. Any (reliable) determination of this contribution will have an 
immediate impact on our understanding of the theoretical errors. On the one hand, it would 
put on more solid basis our implicit assumption that the leading order solution corresponds 
to a Coulomb-type bound state and, once this is achieved, it would move the error estimates 
of the non-perturbative effects from a qualitative level to a quantitative one, (hopefully) 
bringing them down significantly. On the other hand, one may think of cross-checking our 
result with other determinations. The fact that the difference happens to be relatively tiny 
supports our believe that (perturbative and non-perturbative) higher order effects are indeed 
not very large. Alternatively, one may search for combinations of observables less sensitive 
to long distance physics effects in order to get a more accurate result for the masses. 

Leaving aside theoretical errors, one may expect to bring down the errors associated to 
other parameters of the theory significantly (if one has a large enough set of observables) by 
using global fits. 

Another issue that deserves further consideration is whether it is possible to develop a 
renormalon subtraction scheme completely within dimensional regularization (see the discus- 
sion at the beginning of sec. f|). This would provide a better understanding of the physical 
system and therefore the errors could be estimated in a more reliable way. 

We expect to come back to these issues in the near future. 
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